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Abstract 

The production of photon pairs in e + e~ — > 7/77 processes is studied 
using exact helicity amplitudes at tree level. Total cross sections, including 
initial state radiation effects, are given. They are presented in the case of 
quarks as a function of y cnt in the JADE algorithm. In the case of leptons 
we use cuts similar to the ones employed in a recent L3 analysis. Masses of 
the final state fermions are taken into account when appropriate. The cross 
section e + e~ — > t + t~77 is about 10% larger when the r mass is neglected. 
We obtain, with a different method, results which are in good agreement 
with L3 Montecarlo simulation. 
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Introduction 



Photon emission is of crucial importance in e + e~ collisions at energies close to the 
Z° peak, and a good understanding of its properties is necessary in order to test the 
Standard Model with high accuracy. Most of the radiation is soft and/or collinear. 
The cross section for events with hard, isolated photons is much smaller, being ap- 
proximately reduced by a factor a em for each emitted photon with respect to the cross 
section with only soft radiation. Hard photons are interesting because they probe the 
quark charges and because they could be a signal for new particles or for some kind of 
substructure of known particles |], || . Not surprisingly, it has been shown || that the 
soft approximation is inadequate to describe hard photon radiation, and that a matrix 
element calculation is needed for this purpose. 

The L3 Collaboration Q has recently reported four events in which two hard photons 
are emitted with a total invariant mass close to 60 GeV. From a comparison with the 
YSF3 Montecarlo f|, after full detector simulation, L3 estimates the probability for all 
four events to originate from ordinary QED to be of the order of 10~ 3 . The expected 
rate for the production of a Standard Model Higgs of about 60 GeV, in association 
with leptons, decaying in the two-photon channel, is about four orders of magnitude 
smaller than the measured one. 

The obvious thing to do in order to understand whether these peculiar events repre- 
sent the first manifestation of some new physics is to examine all multi-photon events 
and compare them with the available simulations in search of anomalies. 

The L3 events have a typical y cnt in the 10 -3 range, and events with y cnt as small 
as 10~ 4 would still survive all cuts. This range has only been partially analyzed, with 
a simplified set of cuts, in previous studies ||. 

In two recent papers || 0] we have shown that at the Z° peak, when masses are 
properly taken into account, cross sections involving c or b quarks differ significantly 
from the corresponding predictions obtained when masses are neglected. For final states 
with more than one photon the 6-quark contribution is severely reduced by the small 
electric charge of the b. For c-quarks the mass effect is of the order of several percent. 
However, in order to test this result, c-quarks tagging is required, and this substantially 
decreases an already small rate and increases the sources of uncertainty. These results 
suggest that a calculation for the more easily accessible reaction e + e~ — > r + r~77, in 
which the r mass is not neglected, could give significantly smaller predictions than 
those available so far. 

The main purpose of this paper is to give a realistic estimate of the total cross section 
and of some typical distribution for e + e~ — > t + t~77, taking into account initial state 
radiation and the r mass, with a set of cuts which mimic as closely as possible the 
ones used in the L3 experiment. At the same time we study two-photon production 
in association with jets, for which no anomalous event has been reported. We also 
compute the cross section and some distributions for e + e~ — > fi + fi~ , ~y'~y in order to be 
able to compare our calculations with L3 results and to determine if our approach is 
adequate. 



Calculation 



We consider processes without electrons in the final state, which proceed through e + e~ 
annihilation into a photon or a Z°. Both contributions have been retained in all am- 
plitudes. In our calculation hard photon radiation is emitted only from the final state 
fermions, while the radiation from the initial electron legs is described using the struc- 
ture function approach ||. This separation of final and initial state radiation makes 
sense since the two terms are separately gauge invariant and we are only interested 
in photons at relatively large angle with respect to the beam. We have computed all 
matrix elements both in the formalism of |9|, |H| and in that of ||11|| . The analytic 
expressions we have used, in the formalism of [^|, [TIJ, can be easily extracted from 
Appendix B of 0. 

The amplitudes have been checked for gauge invariance. We have used M z = 91.1 
GeV, Y z = 2.5 GeV, sm 2 (8 w ) = .23, a em = 1/128 for the vertices connected to the 
Z°,7 propagator and a em = 1/137 everywhere else, m T = 1.78 GeV and m b = 5. GeV 
in the numerical part of our work. 

In what follows we neglect all hadronization effects, and apply cuts at the partonic 
level. 



Results 

The total cross sections for e + e~ — > and e + e~ — > r + r _ 77 at y/s = M z = 91.1 

GeV are given in table I. For comparison we also give the cross section which results 
neglecting the r mass. Only photons satisfying m 77 > 2.5 GeV, |p 7 | > 1 GeV and 
|cos# 7 | < 0.9 are accepted. For /i's we require |p M | > 3 GeV, \cos9^\ < 0.81 and 
9 7fl > 5°. For r's and for £ ' s ! the reference massless leptons we compare the r's with, 
the assumed cuts are |cos^| < 0.74 and 6* 7 £ > 15°. 

In fig.l we present the total cross sections for e + e~ — > qcpcy for u, d and 6-quarks 
separately and then for the sum over five flavors as a function of y cn f Only for b- 
quarks the mass is taken into account. For photon pairs in association with jets we 
have required 6 iq > 15° but have not imposed any cut on the angles with respect to 
the beam. All particle pairs must have y > y cu t, where y is the usual JADE variable 

HI- 

In fig. 2 we present the energy distribution of the most energetic photon in e + e~ — > 
^+^77, e + e~ — > r + r^'yy and e + e~ — > i^oll with £ massless. In fig. 3 we give 
the two-photon mass distribution in e + e~ — > r + r~77 and again compare it with the 
results for a massless lepton with the same choice of cuts. Finally in fig.4 we show 
the two-photon invariant mass distribution for e + e~ /i + /i~77, e + e~ — ► t + t~77 and 
e + e~ — ► qq'j'j summed over all flavors. 

In fig. 2 and 4 the results are deliberately presented in a way that allows an easy 
comparison with L3 data. 

In order to discuss the consistency of our results with L3 Montecarlo simulations 
we proceed as follows. From the fraction of muon events in the Montecarlo sample 
which survives all cuts, an efficiency of about 56% can be inferred for the L3 detector. 
The main factors contributing to the efficiency are the geometrical acceptance and 
the effect of the dead region between detector modules in the R-0 plane within the 



fiducial volume |cos# M | < 0.81. The first contribution can be easily extracted from the 
ratio of the integrated Born cross section for for e + e~ — > p + fi~, in the presence of the 
angular and momentum cuts and in their absence, and turns out to be about 74%. 
The additional effect related to inefficiencies within the fiducial volume can therefore 
be estimated at about 76%. 

L3 data include runs from 1990 through July 1992 at various energies around the 
Z° peak. An equivalent peak-luminosity of about 22 pb" 1 can be derived, dividing the 
total sample of 950,000 Z° by the peak-cross section of about 43 nb, including initial 
state radiation. This luminosity and the 76% efficiency in detecting muons can then 
be used to translate the cross sections we present, which already include geometrical 
acceptance, in a prediction for the number of events. As an example we would expect 50 
events with two photons and a yU + p~ pair in the sample. In order to compare this result 
with the L3 Montecarlo simulation it is necessary to correct for the cut in the invariant 
mass of the photon pair which is present in our analysis and is not implemented by L3. 
From the m 77 distribution in fig. 6 of |4| one can read that the first mass bin, which 
exactly corresponds to our excluded region, contains 20% of the Montecarlo events. 
Therefore our prediction of 50 events has to be compared with 45 events predicted 
by L3 in the same mass range. This agreement at the 10% level is, in our opinion, 
quite satisfactory. The integrated cross section for m 77 > 50 GeV is 3.5 x 10~ 2 pb for 
/i + yU~77, corresponding to .6 events. 

The consistency of the results obtained in the two approaches, for photon pairs in 
association with leptons, can be better appreciated looking at the shape of the mass 
distribution in fig.4 and in fig. 6 of Qj. There is also good agreement, for large photon 
energies, for the shape of the energy distribution of the most energetic photon shown 
in fig.2 and in fig.2a of ||. The low energy peak, present in fig. 2a of JO, is strongly 
reduced in our distribution due to the additional minimum invariant mass requirement. 

On one hand we interpret these results as a confirmation of the correctness of L3 
Montecarlo simulation, on the other hand they give us confidence in our predictions. 

From table I we see that the predicted cross section for e + e~ — ► r + r _ 77 is about 
9% smaller than it would be for a massless lepton. Fig.2 and 3 show that the effect of 
the r mass is more pronounced for large photon energies and for large invariant mass 
of the photon pair. As expected, the total cross section is dominated by low energy 
photons and low mass pairs. 

The cross sections for photon pairs in association with jets, with the assumed set of 
cuts and y cut = 5 x 10 -3 , is about 1 pb corresponding to 22 events. Since the mentioned 
cuts are similar to those employed in we expect a number of such events to be 
present in the data sample of each LEP experiment. 

As a byproduct of the calculations presented here, we have verified that the ratio of 
distributions for massive particles to those for massless ones is not modified when initial 
state radiation is included. Therefore, also in this more realistic case, all conclusions 
drawn in || [7| remain valid. Obviously all distributions have to be scaled down by 
approximately 30%. 



Conclusions 

We have produced the total cross section and some typical distributions for e + e~ — > 
t + t~77, taking into account initial state radiation and the r mass, with a set of cuts 
which mimic as closely as possible the ones used in the L3 experiment. The total 
cross section is 9% smaller than it would be for a massless lepton and the discrepancy 
increases for larger photon energies and two-photon invariant masses. Photon pair 
production in association with jets has also been computed and we believe it would be 
important to look for such events. We have studied e + e~ — > /i + /i~77 obtaining a total 
rate and distributions in good agreement with L3 results. 
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Table Captions 



table I. Cross sections, cuts and parameters for e e — > 11 77 with I 
where £q is a massless lepton. Errors are as calculated by VEGAS |L3 



/i,r, £ , 



Figure Captions 

fig.l. Cross sections for e + e~ — > uu-yy, e + e~ — > (ic/77, e + e~ — > 6677 and for the sum 
over five flavors as a function of y cut at yfs = 91.1 GeV. The cuts used for photon 
pairs in association with jets are y cut > 5 x 1CT 3 , 9 7q > 15° and |p 7 | > 1 GeV. 

fig. 2. Energy distribution of the most energetic photon for e + e~ — > /i + /i~77, e + e _ — > 
£q£q 77 and for e + e~ — > t + t~77, where £0 is massless, at yfs = 91.1 GeV. Cuts 
as in table I. 

fig. 3. Two-photon invariant mass distributions for e + e~ — > ^o~^o 77 an d for e + e~ — » 
r + r _ 77, where £ is massless, at s/s = 91.1 GeV. Cuts as in table I. 

fig. 4. Two-photon invariant mass distributions for e + e~ — > yU + /i~77, e + e~ — > t + t~jj 
and for e + e~ — > qqy'J, with y cnt = 5 x 10~ 3 and summed over five flavors, at 
y/s = 91.1 GeV. Cuts as in table I for leptons. The cuts used for photon pairs in 
association with jets are y cut > 5 x 1CT 3 , 6 iq > 15° and \p 7 \ > 1 GeV. 



cr(pb) 


e + e — > 77 


e + e — > t + t 77 


e+e -> £+£ 77 


(4.093 ±0.016)* 


(1.4746 ±0.0038)* 


(1.6169 ±0.0052)* 


(3.025 ± 0.015)+ 


(1.0854 ±0.0035) t 


(1.1913 ±0.0048)+ 


\p e \ > 3 GeV 
\cos9 e \ < 0.81 e yt > 5° 


|cos0 € | < 0.74 


9 je > 15° 


m 77 > 2.5 GeV 


|p 7 | > 1 GeV 


Icose^ < 0.9 


v /i= M z = 91.1 GeV 


T z = 2.5 GeV 


sin 2 9 W = 0.23 


* Without ISR; t With ISR. 


Table I 





